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Abstract 

Labour productivity distribution (dispersion) is studied both theoretically and empiri- 
cally. Superstatistics is presented as a natural theoretical framework for productivity. The 
demand index k is proposed within this framework as a new business index. Japanese pro- 
ductivity data covering small-to-medium to large firms from 1996 to 2006 is analyzed and 
the power-law for both firms and workers is established. The demand index k is evaluated in 
the manufacturing sector. A new discovery is reported for the nonmanufacturing (service) 
sector, which calls for expansion of the superstatistics framework to negative temperature 
range. 
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I. Introduction 



Standard equilibrium theory in economics implies that the labour productivity is equal 
among firms and sectors. This was shown to be wrong by a detailed study of the real data 
(Aoyama et al. 2008b[). Facing this situation, one may argue that (i) since the real data is 



slightly different from ideal case that the economic theory deals with and is contaminated 
with various inaccuracies and errors, so slight deviation from the theoretical prediction is 
unavoidable and is even expected, and besides, (ii) since the equilibrium theory is self- 
consistent, reasonable and convincing it must be true, These claims are not valid, as (i) the 
distribution of productivity is wide-spread; it is not even a normal distribution or log-normal 
distribution as expected from contamination argument, but does obey Pareto law (power 
law) (Pareto 1896) for large productivity, that is, the distribution has the distinct char- 
acteristics of fat tails, and (ii) there may exist other theories that are far more convincing 
and the validity of the theory can be judged only by facing the true nature of the subject. 
Indeed, physics, or any other discipline of exact science managed to develop to the current 
status just by following (ii): No matter how the pre-Copernicus theory is reasonable, beau- 
tiful and convincing, earth moves; no matter how the idea of absolute time in the Newtonian 
mechanics, which by the way underlies the current equilibrium theory of economics, seems 
unavoidable, Einstein's relativity theory describes the true nature of time and space. These 
and other numerous historical examples in exact science teaches us that we need to face the 
phenomena seriously and has to construct the theory that meets its demand. Simply put, 
we need to take scientific approach. 



Such was the thought behind the study of the productivity by Aoyama et al. (2008b) 



who proposed the superstatistics theory in statistical physics as the theoretical framework 
for the productivity. 

In this paper, we further advance the superstatistics theory of productivity by examining 
the whole spectrum of firms in Japan, while in the previous work of Aoyama et al. (|2008b) 



and Aoyama et al. (2008a I the data was limited to listed firms. Furthermore, we anal- 



yse the manufacturing and nonmanufacturing sectors separately. In Section II (and in the 
Appendix), we present the superstatistics framework for the productivity for completeness. 
Then in Section III, after explaining the nature of the database and the method of analysis, 
we present the results for the manufacturing sector and the nonmanufacturing sector sep- 
arately. We also study the distribution of the productivity of business sectors. Section IV 
contains some conclusions and discussions on the necessities of extending the superstatistics 
framework. 



II. Superstatistics theory of productivity 

We first review the superstatistics theoretical framework of productivity proposed by | Aoyama 



eTaT1 ( p008bl a]) in a concise manner. The data analysis and the discussion of the evaluation 



of the property of the aggregate demand is done in later sections using this framework. 



A. Statistics 



Yoshikawa and Aoki (Yoshikawa 2003 Aoki and Yoshikawa 20071 proposed an equilibrium 



theory of productivity distribution several years ago. Its essence is the equilibrium theory 
statistical physics, where the most common distribution is realized under given constraint. 
Its beauty lies in the fact that it does not depend the details of the individual properties and 
interactions among constituents (firms in economics and atoms and molecules in statistical 
physics). Let us first review it very briefly 
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We label firms by an index k = 1,2,... ,K, where K is the total number of firms, the 
number of workers at the fc-th firm by n^, the productivity of the fc-th firm by all for a 
given year. There are two constraints on these quantities. 

(i) Total number of workers N: 

K 

J2n k = N. (1) 

k=l 

(ii) The aggregate demand D: The sum of firm's production is the total production, 

which is equal to the aggregate demand D; 

K 

n kc k = £>. (2) 

k=l 

Implicit here is that we are dealing with the mean labour productivity 

c:=|, (3) 

where Y is the value added and L is the labour (in number of workers). Although this 
is different from the marginal labour productivity c ( mar s inal ) := OY/dL relevant in the 
standard equilibrium theory, this difference is irrelevant as we will elaborate later. 

By using the standard proposition that the distribution that maximizes the probability 
under these constraint is realized in nature, which is equivalent to the entropy-maximization, 
the Boltzmann law is obtained. This states that the probability pk of the worker's produc- 
tivity being equal to cj, is the following. 

Pk ■= ^ = — e-^* (4) 



where Z(f3) is the usual partition function: 



K 



Z(/3):=5>-^. (5) 



fe=i 

This guarantees the normalization of the probability pk ; 

X> = 1. (6) 

fc=i 

The parameter is inverse-temperature determined by the mean demand D as follows: 

In our database, we have nearly half a million firms and about 10 million workers (see 
Figj3|. Therefore, it is most appropriate to use the continuous notation, in which the 
probability distribution function (pdf) of the firm's productivity is denoted by p^ F \c), and 
the pdf of the worker's productivity by p( w \c). From Eq.(k|, they satisfy the following: 



P (W) ( C ) = ^ye-^ F >( C ), (8) 

where the partition function is 

Z[p) := / e- pc p^\c)dc. (9) 
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B. Superstatististics 

Although the theoretical prediction ^ is both elegant and powerful, it is quite limited in 
the sense that it is realized in a stationary environment. Namely, the demand (and thus 
the temperature) has to be constant. However, the demand is rarely constant. Rather it 
is one of the most quickly changing parameter (Yoshikawa 20031. Therefore, we need to 
expand the horizon of the theory to meet the changing environment. Just such a theory, 
named superstatistics (statistics of statistics) has been proposed recently by Beck and Cohen 
(20031 in the context of statistical physics. In this theory, the system goes through changing 
external influences, but is in equilibrium described by Boltzmann distribution ^ within 
certain limited scale in time and/or space. Therefore, the whole system can be described by 
an average over the Boltzmann factors, with the weight given by the relative scales (in time 
and space) of the temperature which the system experiences. 

This superstatistics was successfully applied to various systems ( Beck 2005 2008 ) . Most 



analogous to our economic system of firms and workers maybe the Brownian motion of a 



particle going thorough changing temperature and viscosity ( Ausloos and Lambiotte 2006 



Luczka and Zaborek| |2004| ). Our workers are the particles in Brownian motion: They move 



from firms to firms, which keeps trying to meet ever-changing demand by employing and 
dismissing workers. Therefore, the superstatistics is the right framework to deal with the 
distribution of the workers. The weighted average over the temperature now replaces the 
Boltzmann factor e~@ c : 



B{c) 



e- 0c f/3(j3)d(3. 



(10) 



In this equation the changing environment is represented by the weight factor //?(/?), which 
is, in turn, is a function of the mean demand D by Eq.(|7]). The pdf of worker's productivity 
(fsl) is now modified to; 



P (W) (c)=^ (F) (c)S(c). 
The new partition function Zb in the above is given by, 



(11) 



P 



(F \c)B(c)dc. 



(12) 



Let us study what the superstatistics theory tell us for the high productivity region. We 
concentrate on this region because, as we will see by the data analysis in the next section, 
both the firm's productivity and worker's productivity obeys the Pareto's law (power law): 



p (w) (c) oc a'™' 1 . 



(13) 
(14) 



This feature brings advantage to the high-productivity study because of the following rea- 
sons: 

(i) This feature is quite evident in the data and the Pareto indices /xf,w can be estimated 
reliably. In comparison, medium to low range sometimes shows two-peak structure, 
which makes it difficult to extract notable, representative features. (Elsewhere in this 
volume, Souma et al. (20091 elaborates on this point.) 

(ii) As was proven by Aoyama et al. (2008b I, if the Pareto law holds for the "mean" 



productivity c, the same law with the same value of the Pareto index holds for the 
marginal productivity c ( mar s mal ) under a wide assumption. 
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Let us now study the behaviour of Eq.(10l for large c. This integration is dominated 



by the small (3 region. Thus the behaviour of the pdf fp((3) for small (3 is critical. Let us 
assume the following in this range: 



fM « /r 7 ( 7 <i), 



(15) 



where the constraint for the parameter 7 comes from the convergence of the integration in 



Eq.(10l. This leads to the following for large c: 



B(c) ocL(l - 7 )c 



7-1 



(16) 



Substituting this and the Pareto laws Eqs. (13 1 and (14 1 into Eq. (Ill, we obtain the 
following: 

/Xw = /iF — 7 + 1- ( 17 ) 



We note here that because of the constraint 7 < 1 , this leads to the inequality 

Mw > Mf- 



(18) 



This becomes a critical test of this superstatistics theory of productivity, which we come 
back to in the following section. 



The above derivation of Eq.(17l proves, in effect, that the Pareto law for firms and that 



for workers are compatible only if the temperature distribution obeys Eq. (15), as no other 



behaviour could result in the power law ( 16 1. In this sense, we see that empirical observation 



leads to Eq. (fbeta) . This in turn leads to empirical laws to the distribution and fluctuation 
of the demand through Eq. d7l: In this manner, the parameter 7 in the distribution of (3 
is related to a parameter in distribution of D 7 which we denote by fo(D) by the following 
relation: 

fp(f3)d{3 = f D (D)dD, (19) 
Mathematical relation between (3 and D is studied in detail in the appendix. Using the 



result (43 1 and Eq.(15l, we find that 



f D (D) cx «c) - D) 



-s 



with 



_ J S - 1 for 2 < fj, F ; 

7 ~~ \(fi F - 1)(6- 1) forl</x F <2. 



(20) 



(21) 



Note that /d(-D) has an upper limit (c)o: From Eq.([8|, it is evident that as the temperature 
1/(3 goes up, workers move to firms with higher productivity. As the temperature becomes 
infinity (3 = 0, all the firms has the same number of workers. Thus the total demand is 
limited by the values achieved at this point, where D = (c)q. 

Combining Eqs. ( 17 1 and ( 21 ) , we reach the following relation between the Pareto indices: 



-6+1 



for 2 < /ip 



(^ F - l)(-<5 + 1) + A*f for 1 < hf < 2. 



(22) 



This relation between nw and fip is illustrated in Fig[T] The range of the parameter S is 
—00 < 6 < 1 from the normalizability of the distribution of /d(I?). The upper limit may 
also be obtained from the constraint 7 < 1 and Eq. pTj ) . Because of this, Eq. ( 22 ) predicts 
that /xw is larger than [If- Also, Eq.(22l has a fixed point at (/xwjMf) = (1, 1); the line 
defined by Eq. ( 22 1 always passes through this point irrespective of the value of 8. The Pareto 
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Figure 1 : Illustration of the relation between //\y and /iF ( 22 ) . The solid line is 
the relation (22), and the filled circle is the data. 




index for firms is smaller than that for workers, but it cannot be less than one, because of 
the existence of this fixed point. 

This way, the parameter S calculated from fi-p and /z\y represents the behaviour of the 
demand close to its upper limit. As a parameter with the same function, we propose the 
following parameter, which we call Demand Index: 

n := ^ . (23) 

This parameter is a monotonically increasing function of S and ranges from to 1. The 
limited range of k makes easy to handle and plot. If k is close to one, the demand fluctuated 
to the high region significantly; if it is equal to zero, the demand does not go very high (it 
could be dumping faster than any power law toward the upper limit). 

More generally, a function (A — 1) /(A — 5) with A > 1 has the required property, but by 
choosing A = 2 we obtain dn/d5\s=i = 1, so that k ~ 6 for 8 — ► 1. This proximity of k and 
6 is desirable to some extent as the data often shows 6 in the range 0.5 to 1. 

In summary, the superstatistics framework predicts that the Pareto indices /i\y and fip 
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Figure 3: Total number of firms (left) and workers (right) in our database, 
both for manufacturing sector (dashed lines), nonmanufacturing sec- 
tor (dotted lines) and both sectors (solid lines). 



2006 



determines the Demand Index k as follows: 

1 



/Av - Mf + 1 



Hw - 1 

This relation is illustrated in Fig(2j 



for 2 < fi F ; 



for 1 < fi F < 2. 



(24) 



III. Empirical Facts 

A. Database and the Analysis 
In calculating the productivity c by Eq. (|3| from data, we calculate the value added Y by the 



method put forward by the Bank of Japan ( |Souma et al. 20091, which is the most common 
method used in Japan. As for the number of workers L, we use the average of the value of 
that year and that of the past year, as each are defined to be the value at the end of the 
year. 

For comprehensive, high-accuracy study, we made a database from two sources: 

Nikkei-NEEDS Nikkei Economic Electronic Databank System (NEEDS) database is a 
commercial product available from Nikkei Media Marketing, Inc.| p008 1 and contains 
financial data of all the listed firms in Japan. This is a well-established and represen- 
tative database, widely used for various purposes from research to practical business 
applications. We have extracted data from their 2007 CD-ROM version, which contains 
two to three thousand firms and five to six million workers. 



CRD Credit Risk Database (Credit Rise Database Association 20081 is the first and only 



database for small-to-medium firms in Japan. It started collecting data from both 



banks and credit guarantee corporations ( Credit Guarantee Corporations 2008 1 since 
2001. The latter, however, does not contain enough database entry necessary for the 
Bank-of-Japan method and thus are omitted from our database. 

By combining these two database and removing any overlap, we have obtained a unified 
database that covers a wide range of firms in JapanQ The total numbers of the firms and 
workers covered in this unified database is plotted in Fig|3] 



1 The same database was used for productivity analysis by Ikeda and Souma (20081 
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Unfortunately, CRD data covers only from 1996 to 2006, and does no go far in the 
past, unlike Nikkei-NEEDS. This limits our unified database to the same short period. We, 
however, consider that it is important to include small-to- medium firms to our analysis so 
that we have a view of the whole spectrum of firms and workers in Japan. Thus we have 
decided to create and analyse this database. In coming years, we trust that CRD will keep 
collecting data. So, we may extend our analysis to future, if not past, by continuing and 
expanding what we have started here in this paper. 

We need a model distribution to fit the data and extract the value of the Pareto indices. 
It has to have several properties: (i) It has be defined in the whole range, c € [0,oo]; (ii) it 



has to have power law (13 1 and (14 1 for large c; and (iii) it has to be able to describe the 
data to some accuracy over the whole region. The power law manifest itself as the straight 
line in the log-log rank-size plots of the data (see Fig 4[^) Since its gradient is the Pareto 
index, its value is estimated fitting the straight section of the rank-size plot with a straight 
line. Although this can be done easily and is intuitive, it has several pitfalls: Often, the 
definition of "the straight section" is ad-hoc. Slight change of it can bring nonnegligible 
change in the value of the Pareto index. Even if a good one can be found for a particular 
year, it might not work for other years of the same database, which make comparison of 
different years meaningless. 



The "Generalized Beta Distribution of the Second Kind" (GB2) (Kleiber and Kotz 20041 



satisfies the property (i)-(iii) and yet manageable. It is defined by the following pdf; 



P 



(GB2J 



9 



1 / c 



1 



where the four parameters satisfy constraints fx, v, q, c\ > 0. Since for large c; 



P 



(GB2) 



(C) ~ - 



Cl 



(25) 



(26) 



the parameter \i is the Pareto index. Incidentally, its cumulative distribution functions (cdf) 
is the following: 



P 



(GB2) 



(c) 



P 



(GB2) 



(c) dc. 



B(z,n/q, v/q) 
B((J./q,v/q) 



1 + 











(27) 



where B(z, s, t) is the incomplete Beta function with 5(1, s, t) — B(s, t). (Detailed study of 
small-to-medium productivity was done by Souma et al. (2009 1 using this GB2 distribution.) 



B. Manufacturing firms 

The rank-size plots of the productivity of the manufacturing sector in 2004 is given in Fig]4] 
by dots, together with the best-fit cdf obtained by the maximum likelihood method. In these 
log-log plots, we see that the actual distributions of the data are close to straight lines for 
large c, which implies that it obeys the power law (i.e., the Pareto law) as we have discussed 
above. The best-fit cdf (dashed lines) indeed represents the data to good accuracy. The 
situation is quite similar in all other years. 

The values of the Pareto indices /ip and /i\y thus obtained are plotted in Fig(5] Sub- 
stituting these values to Eq. ( 24 1 , we have obtained the value of the demand index n joined 



by solid lines in Fig|6] We see here that the demand is slowing rising during this period, 



2 Firms with extremely high value of productivity are removed from this plot, as they often report one worker, 
which, in view of their huge income, cannot be a good representation of their manpower. 
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Figure 4: The rank-size plot of the productivity data (dots) and the best-fit cdf 
(dashed line) for firms (left) and workers (right) in 2004. 




Figure 5: The Pareto indices z/p and ixw for the manufacturing sector. 
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Figure 6: The demand index k for the manufacturing sector (solid lines) and 
the Nikkei Business Index (dashed lines). 
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Figure 7: The Pareto indices /xp and //\y for the nonmanufacturing sector. 



which is in agreement with general observations in Japan. Plotted in Fig|6]with dashed lines 
is the Nikkei Business Index (NBI), which is a major business index in Japan (Nikkei Net 



Interactive 20081. We observe here that their correlation is good to some extent, which is 



consistent with the fact that our demand index n provides a measure of demand. 



C. Nonmanufacturing (service) firms 

In the nonmanufacturing sector, the same analysis leads to the result plotted in Fig(7] It 
is quite notable that it is completely different from the manufacturing sector: The Pareto 
index /ip is larger than /i-yv- Since the larger Pareto index means that the pdf is damped 
highly for large c, this means that the higher productivity firms, more workers are employed. 
This is not allowed under the ordinary Boltzmann distribution (JsJ) due to the Boltzmann 
factor e~ l3c . It is not allowed in the superstatistics either, since it is an weighted average 
over the Boltzmann distribution. Therefore, this behaviour of the nonmanufacturing sector 
calls for extension of the theoretical framework. 
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Figure 8: Productivity distribution of 26 sectors from 1996 to 2006. The dashed 
line is a power-law distribution with the Pareto index fj, = 1.6. 



IV. Productivity of Sectors 

Productivity of sectors are of interest. Our database, following Nikkei NEEDS, contains 
26 sectors, among which 15 are manufacturing sectors and 12 nonmanufacturing sectors. 
Their productivity distributions from 1996 to 2006 is plotted in FigjS] Evident in this plot 
is that the productivity distribution obeys the power-law with the Pareto index /i ~ 1.6 
every year. Since the number of data is limited, unlike the firms and workers, fitting with 
GB2 distribution and estimating the value of fi is not very illuminating. In other words, the 
obtained values of fi would suffer from large statistical errors. It is more so if manufacturing 
sectors and nonmanufacturing sectors are studied separately. 

The notable feature of the sector distribution is the (i) it is approximately Pareto and 
(ii) the values of the Pareto index are certainly lower than that of firms every year. This is 
in accordance with the idea of applying superstatistics framework to firms and sectors, in 
contrast to workers and firms ( Aoyama et al.[ [2008b[ ): We may now think of a firm (instead 



of a worker) choosing a sector (instead of a firm) for its business activity. Applying the 
superstatistics to them, we find again that the Pareto index of the firms are higher than 
that of sectors. This is what we observe in FigM 



V. Conclusion and Discussions 



We have studied the superstatistics theory of productivity and have proposed the demand 
index k, which determines the relation between the Pareto indices of the productivity dis- 
tributions of firms and workers. Analysis of the whole spectrum, from small to large, firms 
in Japan from 1996 to 2006 is carried out and manufacturing sector was studied within the 
superstatistics framework. 

One might argue that what we have been observing is a temporal situation and even- 
tually, the effect of shocks, including 1999 bubble collapse, would be averaged out and the 
productivity starts converging to a unique value, as predicted by the orthodox equilibrium 
theory. In order to study this point, we have examined the change of the productivity of 
sectors over the years. The result is plotted in Fig(9j Each line represents a sector, solid 
lines for manufacturing and dashed lines for nonmanufacturing. In this Figure, we see that 
they are far from arriving at a unique value. Rather, they keeps fluctuating. Sometimes, the 
difference of productivity between the sectors widens. This sort of behaviour is typical in 
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Figure 9: Evolution of productivity of each sector, sectors, (a) is the actual 
value, while (b) is the value normalized at 1996. Solid lines are of 
manufacturing sectors and dashed lines are of nonmanufacturing sec- 
tors. 



physical systems. The distribution is maintained, while viewed in detail, each atoms keeps 
changing its energy and momentum. Such is the nature of physical equilibrium and so is 
the economic system. 

We stress here that the distribution ( 20 1 of the demand D is unique in the superstatistics 
framework: It has to obey this distribution if Pareto law hold for both firms and workers 
and Eq. ( 18 ) is satisfied. Just what kind of dynamics and economic principles underlies 
the demand distribution (20 1 is unknown. It would be quite interesting and important to 
construct a model which leads to this kind of behaviour. It is quite possible that such a model 
could be one of the building blocks useful and necessary for reconstructing macro-economics. 

On the other hand, nonmanufacturing (service) sector showed peculiar characteristics 
that has never been seen before. The Pareto index for the workers were larger than that for 
firms. 

In the ordinary Boltzmann distribution ( 10 ) , the positivity of the temperature guarantees 
that higher the productivity less workers are employed. Since superstatistics is the weighted 
average of the Boltzmann distribution, no matter what the weight function fp(f3) is, the 
higher productively means less workers, which is the reason for fi-p < /i-yv- The fact that the 
nonmanufacturing sector violates this constraint means that the weighted average over the 
negative temperature is required. 

Negative temperature is possible for a physical system in nonequilibrium. One such an 
example is a laser, where many atoms or electrons are in a excited state before the emission. 
For the current case, there is at least one economic reason why it is required; excess of the 
demand. There are certain limits to the productivity of a given firm due to many constraints 
it faces. But by hiring more people while maintaining the same structure, firms can increase 
its add value, thereby meeting the increasing demand. Just such cold be happening in the 
nonmanufacturing industry. 

On the other hand, superstatistics of the negative temperature has not been developed 
yet: While it is easy to expand the weighted average over the negative temperature, its full 
consequences are not clear at this stage. It would therefore be quite interesting to develop 
this theory and use it to deal with the nonmanufacturing sector. 
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Appendix: Temperature and the Demand 

We first note the following three basic properties (i)-(iii). 



(i) The temperature, T = 1/(3 is a monotonically increasing function of the aggregate 
demand, D. We can prove it using Eq.([7]) as follows: 

%—h%-^ m - f( ^'-^°- (28) 

where (c n )p is the n-th moment of productivity defined as follows: 

1 r°° 

{cn) ^mJ cV^e-^da (29) 

Note that (c)p = D. This is a natural result. As the aggregate demand D rises, workers 
move to firms with higher productivity. It corresponds to the higher temperature due 
to the weight factor e _/3c . 

(ii) For (/3->oo), 

D -> 0. (30) 

This is evident from the fact that in the same limit the integration in Eq.Q is domi- 
nated by c ~ due to the factor e~ f3c . 

(iii) For T — * oo (/3™*0), 

/ cp^(c)dc(=(c) ). (31) 



This can be established based on the property (i) because D = (c)^ — > (c)o as (3 — ► 
and Z(0) = 1. 

Let us now study the small (3 (high temperature) properties. One possible approximation 
for Eq.(|9]) is obtained by expanding the factor e~^ c = 1 — (3c + ■ ■ ■ and carrying out the 
c-integration in each term. This leads to the following: 

Z((3) = ^ p^(c) (l-/3c+ ^(/3c) 2 + ..}jdc 

= l-(c) /3+i(c 2 )o/? 2 + ..., (32) 
where we have used the normalization condition, 

p ( - F \c)dc=l. (33) 



The result (32 1 is, however, valid only for pp > 2 since (c 2 )o is infinite for /xp < 2, which is 
true as we have seen. 

The correct expansion for 1 < fip < 2 is done in the following way. We first separate out 
the first two terms in the expansion of the factor e"^ c ; 

/>oo 

Z((3) = / (c) (1 - (3c + (e-^ -\ + (3c))dc 

Jo 

= l-(c) (3 + Z 2 (f3), (34) 
Z2((3)^J o p {F \c)g{c)dc= (--Pl F) (c))g(c)dc 

'P?\c)^dc, (35) 
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where g(c) = e @ c — 1 + /3c is a monotonically increasing function of c with 

5 (0) = </(0) = 0. (36) 



The c-integration in Eq.(35 1 is dominated by the asymptotic region of c for small /3. There- 
fore, the leading term in Z 2 (P) is evaluated by substituting the asymptotic expression of 

p <r >W; 



into Eq.(35). We thus arrive at the following: 



Co J dc 

= AlF r(- AlF )(c /3)' IF + --- . (38) 



The case /i F = 2 can be obtained by taking the limit /i F — > 2+ in the following expansion 
valid for 2 < /i F < 3: 

Z(/3) = 1 - (c) (3 + i((c 2 ) - (c) 2 )/3 2 + / x F r(-/x F )(c /3)^ + . . . . (39) 

which can be obtained in the manner similar to the above. The third term is finite for 
fiF > 2, but diverges as // F — * 2+ as 

< C 2 )o-^V (40) 

This cancels the divergence of the fourth term in the same limit and the remaining leading 
term is as follows: 

Z(/3) = 1 - (c) /3 ~ (c /3) 2 log(c /3) + . . . . (41) 
In summary, the partition function behaves as follows: 

ri-(c) /3+i(c 2 ) /3 2 + ... for2< MF ; 
Z(P) =\l (c) P - (c /3) 2 log(c /3) + . . . for /i F = 2; (42) 
[l-(c)o / 3 + /x F r(-^ P )(co/3)" p + ... for 1< fi F < 2. 

Substituting the above in Eq.Q, we obtain the following: 

f (c) - ((c 2 ) -(c)g) /?+.... for2< MF ; 
D= i ( c > + 2c 2 /31og(co/3) + ... for MF = 2; (43) 

[ (c) - /i F 2 r(- A/F )c^^^ 1 + . . . for 1< M F < 2. 
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